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D . Abstract 

' We study the problem of preemptive scheduling of n equal-length jobs with given release 

, times on m identical parallel machines. The objective is to minimize the average flow time. 

■ Recently, Brucker and Kravchenko proved that the optimal schedule can be computed in 
polynomial time by solving a linear program with 0{n^) variables and constraints, followed by 

^/J • some substantial post-processing (where n is the number of jobs.) In this note we describe a 

, simple linear program with only 0(rnn) variables and constraints. Our linear program produces 

^ ' directly the optimal schedule and does not require any post-processing. 

1 Introduction 

> 

I In the scheduling problem we study the input instance consists of n jobs with given release times, 

^\ • where all jobs have the same processing time p. The objective is to compute a preemptive schedule 

of those jobs on m machines that minimizes the average flow time or, equivalently, the sum of 
completion times, X^Cj. In the standard scheduling notation, the problem can be described as 

■ P|rj, pmtn,pj = p\ ^Cj. Herrbach and Leung j3| showed that, for m = 2, the optimal schedule 
can be computed in time 0(?7-logn). Du, Leung and Young |5] proved that the generalization of 
this problem where processing times are arbitrary is binary NP-hard. We summarize these results 
in Tabled 

Very recently, Brucker and Kravchenko [J gave a polynomial-time algorithm for any number m 
of machines. Their algorithm consists of two stages: first, they solve a complex linear program with 
' 0(?i^) variables and constraints, which is followed by a post-processing stage where they construct 

an optimal schedule from the optimal solution of this linear program. 

In this note, we show that there is always an optimal schedule in a particular form, which we 
call normal. We then give a simple linear program of size 0{mn), which directly defines an optimal 
normal schedule. As a side-product, we show that there are optimal schedules with only 0{mn) 
preemptions, improving the 0{n^) bound on the number of preemptions from 
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Problem 



Complexity 



P2|rj ; pmtn;pj = p\Yl Cj 
P\rj;pmtii;pj =p\YlCj 
P|r'j;pmtn lYl^j 
P\ pmtn; pj = p\J2Cj 

Ph; Pj=p\ECj 

P2\rj;pnitn;pj = p\J2 '^j^j 
P|rj; pmtn;pj = p\Yl, "^jCj 



solvable in time 0(n log n) [3] 

solvable in polynomial time improved in this paper 

binary NP-complete [2] 

solvable by the greedy algorithm (trivial) 

solvable by the greedy algorithm (trivial) 

open 

unary NP-complete 



Table 1: Complexity of related scheduling problems. P2 stands for the 2-machine problem. In 
problems with the objective function ^ WjCj, each job j comes with a weight Wj, and the goal is 
to minimize the weighted sum of all completion times. 

2 Structural Properties 

Basic definitions. Throughout the paper, n and m denote, respectively, the number of jobs 
and the number of machines. The jobs are numbered 1,2, ... ,n and the machines are numbered 
1,2, ... ,m. All jobs have the same length p. For each job j, Vj is the release time of j, where, 
without loss of generality, we assume that = ri < ... < r„. 

We define a schedule to be a function which, for any time t, determines the set X[t) of jobs 
that are running at time t. This set X{t) is called the profile at time t. Let X~^{j) denote the 
set of times when j is executed, that is X~^{j) = {t : j £ X(t)}. In addition we require that X 
satisfies the following conditions: 

(si) At most m jobs are executed at any time, that is |'V(t)| < m for all times t. 

(s2) No job is executed before its release time, that is, for each job j, if t < rj then j ^ X{t). 

(s3) Each job runs in a finite number of time intervals. More specifically, for each job j, X~^{j) 
is a finite union of intervals of type [s, t). 

(s4) Each job is executed for time p, that is |A'~^(j)| = p. 

It is not difficult to see that condition (s3) can be relaxed to allow jobs to be executed in infinitely 
(but countably) many intervals, without changing the value of the objective function. 

By Cj = sup X^^{j) we denote the completion time of a job j. In this paper, we are interested 
in computing a schedule that minimizes the objective function Cj. 

Note that, since we are dealing with preemptive schedules, it does not matter to which specific 
machines the jobs in X(t) are assigned to. When such an assignment is needed, we will use the 
convention that the jobs are assigned to machines in the increasing order of indices (or, equivalently, 
release times): the job with minimum index is assigned to machine 1, the second smallest job to 
machine 2, etc. 

Some observations. We now show that, for the purpose of minimizing our objective function, 
we can restrict our attention to schedules with some additional properties. 

Call a schedule X left- adjusted if it satisfies the following condition: for any times s < t, where 
|r^(,s)| < m, if j G X{t) and rj < s then j € X{s) as well. Any optimal schedule is left-adjusted. 
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for otherwise, if the above condition is not satisfied, we can move a sufficiently small portion e > 
of j from the last block where it is executed to the interval [s, s + e), obtaining a feasible schedule 
in which the completion time of j is reduced by e and other completion times do not change. Thus 
we only need to be concerned with left-adjusted schedules. 

We say that the completion times are ordered in a schedule A', if Ci < C2 < . . . < C^. Brucker 
and Kravchenko showed that any schedule can be converted into one with ordered completion 
times, without increasing the objective value. 

Irreducible schedules. We say that a schedule X is irreducible if it is left-adjusted and satisfies 
the following condition for any times s < t: 

max{X{s) - X{t)) < mm{X{t) - X{s)). (1) 

In the formula above we use the convention that max(0) = —00 and min(0) = -|-oo, so holds 
whenever X{s) C X{t) or X{t) C X{s). 

For a schedule X and a job j, define the halfway point of j as^ Hj{X) = ^ Jx-^^^j) Then 
let H{X) = {Hi{X),H2{X), . . .,Hn{X)). For two different schedules X, y, we say that H{X) is 
lexicographically smaller than H{y), if Hi{X) < Hi{y) for the smallest i for which Hi{X) 7^ Hi{y). 

Lemma 1 Let X he a schedule and [s, s -\- e), [t, t + e) two time intervals such that s + e <t, for 
some e > 0. Suppose that there are two jobs i < j with rj < s, such that [s, s-|-e) C X~^{j) — X~^{i) 
and [t,t -|- e) C X~^{i) — X~^{j). Let y denote the schedule obtained from X by exchanging jobs 
i,j in intervals [s, s + e), [t, t + e). Then H{y) is lexicographically strictly smaller than H{X). 

Proof: Clearly, Hi{y) = H,{X) - {t - s)e/2, and Hk{y) = Hk{X) for all k < i. This directly 
implies the lemma. □ 

We need to show that there exists an optimal irreducible schedule. This is quite easy to show if 
we put some restrictions on the granularity of the schedules, for example if we assume that schedules 
are constant in each unit interval [t, t -|- 1) for t E N. In that case one can show that after finite 
number of exchanges (as defined in the previous lemma) any optimal schedule can be transformed 
into an irreducible optimal schedule. This applies, in particular, to the case when the processing 
time and all release times are integer ^ Theorem 6]. The proof for arbitrary real numbers is more 
difficult, and is based on the following lemma, whose proof appears in the appendix. 

Lemma 2 There exists an optimal schedule X for which the vector H{X) is lexicographically min- 
imum over all optimal schedules. 

Proof: See Appendix 1X1 □ 

Lemma 3 There exists an optimal schedule that is irreducible. 

Proof: Let X be an optimal schedule which minimizes H{X) among all optimal schedules. 
According to Lemma |^ X is well defined. By optimality, X is left-adjusted. 

We claim that the completion times in X are ordered. Towards contradiction, suppose there 
are jobs i < j with Cj > Cj. Let [t,t -|- e) C [Cj,Ci) be an interval where i is scheduled and 

^This is a standard value in scheduling, even though the factor ^ is irrelevant for this paper. 
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[s,s + e) C [0,Cj) a arbitrary interval, where j but not i is scheduled. This contradicts the 
minimality of H{X) by Lemma H 

We claim that X also satisfies Towards contradiction, suppose it does not. Then there are 
time intervals [s, s + e) and [t, t + e) for s + e < t and jobs i < j such that X schedules j but not i 
in [s,s + e) and schedules i but not j in [t,t + e). Then Lemma ^ applies as before and the proof 
is now complete. □ 

We now give a characterization of irreducible schedules that will play a major role in the 
construction of our linear program. 

For a given job j and a time t we partition X{t) — {j} into jobs released earlier and jobs released 
later than j. Formally, <^<j(t) = {i G X{t) : i < j} and <Y>j(t) = {i G X(t) : i > j}, see figure ^ 
The lemma below provides a characterization of irreducible schedules. 

u s t 



jobs released : , , 
before k ■^<k[u) 



jobs released 
after k 



Figure 1: Structure of any irreducible schedule. 



Lemma 4 Let X he an irreducible schedule. Suppose that we have two times u < t and a job k 
such that rk <u. Then: 

(a) IfkeX{u)- X{t) then \X<k{t)\ < \X<k{u)\. 

(b) Ifke X{t) - Xiu) then \X{u)\ = m, \X^k{t)\ > \X>k{u)\, and \X<k{t)\ < \X<k{u)\. 

(c) Ifke X{u) r\X{t) then \X<^k{t)\ < |'^<fcWI and \X^k{t)\ > \X>k{u)\- 

Proof: (a) If there was a j € X^k{t) ~ ^<k{u), this would imply that ma,x{X{u) — X{t)) > k > 
j > min(^(t) — X{u)), contradicting irreducibility. Thus (a) follows. 

(b) Since k G X{t) — X(u) and < u, the assumption that X is left-adjusted implies that 
[-^(tt)! = m. 

We must have Af>fc(n) C for otherwise, the existence of A; G X{t) — X{u) and an 

/ G Xykiu) — Xykit) would contradict irreducibility. The inequality |^%'>fc(M)| < |<^>fc(t)| follows. 
This, the assumption of the case, and = m imply |^%'<fc(M)| > 

(c) We only prove the first inequality, as the proof for the second one is very similar. Towards 
contradiction, suppose |A'<;fc(ii)| < | and pick any i G X^kit) — X^k{u)- Then ri < rk < u 
and i G X(t) — X(u), and so the assumption that X is left-adjusted implies X(u) = m. This, in 
turn, implies that |r^>fc(n)| > so we can choose j G X^ki^) — X^kit). But this means that 
j > k > i and j G X{u) — X{t), and the existence of such i and j contradicts irreducibility. □ 

3 A Simple Linear Program 

Machine assignment. We now consider the actual job-machine assignment in an irreducible 
schedule X. As explained earlier, at every time t we assign the jobs in X(t) to machines in order. 
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that is job j G ^(t) is assigned to machine 1 + |A'<j(t)|. Lemma |3 imphes that, for any fixed j, 
starting at t = rj the value of |<%'<j(t)| decreases monotonically with t. Therefore, with machine 
assignments taken into account, X will have the structure illustrated in Figure ^ 

Call a schedule X normal if for each job j and each machine q, job j is executed on g in a single 
(possibly empty) interval [Cj^q, Sj^q), and 

(1) Cj^q < Sj+i^q for each machine q and job j < n, and 

(2) Cj^q < Sj^q-i for each machine q > I and job j. 

By the earlier discussion, each irreducible schedule is normal (although the reverse does not 
hold.) An example of a normal (and irreducible) schedule is shown in Figure [2 



1 


2 


3 


4 


6 


8 


9 


10 




2 


3 


4 


5 


7 


9 


10 








3 


4 


5 


6 




8 




10 








4 5 


6 


7 


9 






1 ' 


2 ' 


3 U 1 


6 1 


7 ' 


8 '9 ' 


10 



Figure 2: Example of a normal schedule. The processing time is p = 8. 



Linear program. We are now ready to construct our linear program: 



minimize 



E •=! C,,i (2) 



subject to ~'Sj,m < —fj j = 1, . . . ,n 



2^(Cj,<? - Sj,q) =p 3=1, 



,n 



Sj,q - Cj,q < i = 1, . . . , n, g = 1, . . . , m 

- Sj,g~i ^0 j = I,. . . ,n, q = 2, . . . ,m 

Cj^q - Sj+i^q < i = 1, . . . , n - 1, g = 1, . . . , m 

The correspondence between normal schedules and feasible solutions to this linear program 
should be obvious. For any normal schedule, the start times Sj^q and completion times Cj^q satisfy 
the constraints of And vice versa, for any set of the numbers Sj^q, Cj,q that satisfy the 

constraints of (jSj), we get a normal schedule by scheduling any job j in interval [Sj^q, Cj^q) on each 
machine q. Thus we can identify normal schedules X with feasible solutions of Q- Note, however, 
that in X a job j could complete earlier than Cj^i (this can happen when Cj^i = Sj^i.) Thus the 
only remaining issue is whether the optimal normal schedules correspond to optimal solutions of 

©• 

Theorem 5 The linear program above correctly computes an optimal schedule. More specifically, 
min;f Cj = minx^j Cj^i, where the minima are over normal schedules X, and Cj represents 
the completion time of job j in X. 
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Proof: By the correspondence between normal schedules and feasible solutions of ((21), discussed 
before the theorem, we have Cj < Cj^i for all j, and thus the < inequality is trivial. 

To justify the other inequality, fix an optimal irreducible (and thus also normal) schedule X 
and a job j, and let [Sj^g, Cj,g) be the last (that is, the one with minimum q) non-empty execution 
interval of j. Consider a block [s,t) where t = Cj^g. By Lemma |1] all jobs executed on machines 
1,2, ... ,g — 1 in [s,t) are numbered lower than j. Further, by the ordering of completion times, 
they are not executed after t. Thus they must be completed at t as well. Therefore we can set 
[Sj^h, Cj^h) = [t, t), for all machines 1 < h < q, without violating any inequality. This gives a normal 
schedule in which Cj = Cj^i. □ 

4 Final Remarks 

We proved that the scheduling problem P\rj,pmtii,pj = p\^Cj can be reduced to solving a 
linear program with 0{mn) variables and constraints. This leads to a polynomial time algorithm 
more efficient than the one resulting from The question whether linear programming can be 
avoided, and whether this problem can be solved with a combinatorial, strongly polynomial time 
algorithm (whose number of steps is a polynomial function of only m and n) remains open. Our 
characterizations of optimal schedules could be helpful in designing such an algorithm. 

Brucker and Kravchenko ^ showed that there is an optimal schedule with 0{rfi) preemptions. 
Our proof provides a better, 0{mn) bound on the number of preemptions, since in an irreducible 
schedule each job is preempted at most m — l times. (Of course, we can always assume that m < n.) 
We do not know whether this bound is asymptotically tight. It is thus quite possible that there 
exist optimal schedules in which the number of preemptions is 0(n), independent of m. If this is 
true, this could lead to efficient combinatorial algorithms for this problem whose running time is 
even independent of m, perhaps even as fast as 0{n log n). Such improvement would require a much 
deeper study of the structural properties of optimal schedules. Since we use only the existence of 
normal optimal schedules, rather than irreducible schedules, we feel that the problem has more 
structure to be exploited. 

We implemented the complete algorithm (converting the instance to a linear program and 
solving this linear program). It is accessible at http : //www. Iri . f r/~durr/P_rj_pmtn_pjp_sumCj . 
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A Proof of Lemma [2] 



Lemma 1 There exists an optimal schedule y* that minimizes lexicographically H{y*) among all 
optimal schedules. 

Proof: The major difficulty that we need to overcome is that the set of schedules is not closed 
as a topological space, so there could be a sequence of schedules with decreasing values of H{X) 
whose limit is not a legal schedule. The idea of the proof is to reduce the problem to minimizing 
H(X) over a compact subset of schedules. 

Define a block of a schedule A" to be a maximal time interval [u, t) such that (u, t) does not 
contain any release times and A'(s) is constant for s € [u,t). 

For convenience, let r„_(_i to be any upper bound on the last completion time of any optimal 
schedule, say r„+i = r„ + up. Thus all jobs are executed between ri and r„+i. Each interval 
[rj, Tj+i), for i = 1, . . . , n is called a segment. By condition (s3), each segment is a disjoint union of 
a finite number of blocks of X. Also, for each job j, we have Cj = t for the last non-empty block 
[s,t) whose profile contains j. 

A schedule X is called tidy if all jobs are completed no later than at r^+i and, for any segment 
[rj,rj_|_i), the profiles X{t), for t S [rj,rj_|_i), are lexicographically ordered from left to right. More 
precisely, this means that, for any r^ < s < t < rj+i, we have 

min(;f(s) - ;f(t)) < mm{X{t) - X{s)). 

One useful property of tidy schedules is that its total number of blocks (including the empty ones) 
is = nM, where M = X^JT^q (")• Fi^om now on we identify any tidy schedule X with the vector 
X G M'^ whose 5-th coordinate Xb represents the length of the 6-th block in X. 

In fact, the set T of tidy schedules is a (compact) convex polyhedron in M^, for we can describe 
T with a set of linear inequalities that express the following constraints: 

• Each job j is not executed before r^, 

• Each job j is executed for time p. 

For example, the second constraint can be written as Xf, = p, where the sum is over all blocks 
b whose profile contains j. 

Claim 6 ^ Any schedule X can be transformed into a schedule with ordered completion times, 
without increasing the objective function value. 

Suppose for jobs i < j the completion times in X satisfy Cj > Cj. Then there must be a maximal 
time t, such that in [t, Ci) both jobs are scheduled for an equal amount of time. Exchanging both 
jobs in this interval will reorder their completion times. After repeating this process sufficiently 
many times, eventually all completion times will be ordered. See P for details. 

Claim 7 Let X be a schedule in which completion times are ordered and upper bounded by r„+i. 
Then X can be converted into a tidy schedule X' such that 

(a) Cj < Cj for all j (where Cj and are the completion times of j in X and X' , respectively.) 

(b) H{X') is equal to or lexicographically smaller than H(X). 
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Indeed, suppose that X has two consecutive blocks A = [u, s), B = [s, t) where ri < u < s < t < 
Tj+i, and the profile X{u) of A is larger (lexicographically) than the profile A'(s) of B. Exchange A 
and B, and denote by X' the resulting schedule. Let j = min(Af(s) — X{u)) < m.m(X{u) — X{s)). 
Since Cj > t, all jobs in (n) — X{s) are also completed not earlier than at t. So this exchange does 
not increase any completion times. We have Hj{X') < Hj{X) and Hi{X') = Hi[X) for i < j. Thus 
H[X') is lexicographically smaller than H{X). By repeating this process, we eventually convert X 
into a tidy schedule that satisfies the claim. 

We now continue the proof of the lemma. Fix some optimal schedule X* . Let C| denote the 
completion time of a job j in X* . From Claims IHl and [3 we can assume that X* is tidy and 

< C| < . . . < C*. (For the peace of mind, it is worth noting that Claim [3 implies that X* is 
well defined, for it reduces the problem to minimizing Cj over a compact subset T of M-^.) 

Consider a class Tq C T of tidy schedules X such that each job j in X is completed not later than 
Cj. Since X* G Tq, the set Tq is not empty. Similarly as T, Tq is a (compact) convex polyhedron. 
Indeed, we obtain Tq by using the same constraints as for T and adding the constraints that each 
job j is completed not later than at C*. To express this constraint, if in X* the completion time 
Cj of j is at the end of the a-th block in the segment [ri,rj+i), then for each b > (i — 1)M + a 
such that j is in the profile of the 6-th block, we would have a constraint Xh = 0. Note that these 
constraints do not explicitly force j to end exactly at Cj, but the optimality of X* guarantees that 
it will have to. 

Now we show that there exists a schedule y* € Tq for which H{y*) is lexicographically mini- 
mum. First, as we explained earlier, Tq is a compact convex polyhedron. Let Ti C Tq be the set 
of X for which Hi{X) is minimized. Hi is a continuous quadratic function over Tq, and thus Ti 
is also a non-empty compact set. Continuing this process, we construct sets T2, . . . ,T„, and we 
choose 3^* arbitrarily from T„. □ 
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